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Abstract 

The attractive Casimir force acting on a micrometer-sphere suspended 
in a spherical dip, close to the wall, is discussed. This setup is in prin- 
ciple directly accessible to experiment. The sphere and the substrate are 
assumed to be made of the same perfectly conducting material. 
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1 Introduction 

Experiments aiming at testing the theory of the Gasimir effect (pQ; for recent 
reviews see are more numerous than what one might perhaps think. Let 

us here only highlight some examples, starting with the classic experiment of 
Sparnaay [SJ ^ . This experiment tested the Gasimir force between two parallel 
plates, made of chromium steel, chromium, and aluminium. With the exception 
of aluminium (whose problems most likely were due to impurities), the results 
were in good qualitative agreement with the Lifshitz formula [TJ , calculated from 
the assumption of perfect reflecting boundaries. The experimental technique 
was based upon use of a spring balance (sensitivity about 10"'^ dynes), sensing 
the attractive force. The plates were assigned parallel by visual inspection. 

Another well known classic experiment is the one of Sabisky and Ander- 
son 12], dealing with the properties of liquid helium films adsorbed on cleaved 
surfaces of alkaline-earth fluoride crystals at T = 1.38 K. Film thicknesses mea- 
sured by means of acoustic interferometry, were found to lie between 1 and 25 
nm. At thermal equilibrium the film thickness gets a value that is determined 
thermodynamically, given the Lifshitz formula for the Gasimir force as input. 
The results measured were in very good quantitative agreement with the Lifshitz 
expression. 

The modern series of experiments was initiated with the seminal work of 
Lamoreaux . He used a balance based on a torsion pendulum to measure the 
Gasimir force between a gold coated spherical lens (radius about 12 cm) and a 
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flat plate. The lens was mounted on a piezo stack and the plate on one arm 
of the torsion balance. The Casimir force would result in a torque, which was 
detected via a capacitance measurement. Maximum separation between the two 
surfaces was 12.3 fim. In a later note, Lamoreaux included corrections, such as 
those arising from finite conductivity ^], and with Buttler ^T] he gave recently 
an analysis of thermal noise in torsion pendulums. The Lamoreaux experiment 
gave rise to a surge of activity, both experimentally and theoretically. 

In the most recent years, the atomic force microscope (AFM) used in partic- 
ular by Mohideen et al. has led to the most accurate determination 
of the Casimir force between a micrometer-sized sphere and a plate. By using a 
sphere/plate configuration, one avoids the strict requirement about parallelism 
that is so demanding in the case of parallel plates. The accuracy is now of 
the order of a few per cent; this accuracy being actually under current debate 
mainly because of the temperature corrections. 

We shall not here go into further detail as regards the experimental status. 
A detailed exposition on the experiments up to 2001 is given in the review of 
Bordag et al. referred to earlier ^ , and a detailed survey of the developments 
in the last four years is given in Milton's review [5], Sect. 3.6. We mention, 
though, the impressive plane-plate experiment of Bressi et al. |15| : they were 
able to guarantee a parallehsm of the plates to better than 3 x 10~^ rad. (It 
may even be that this experiment has been the first to measure the temperature 
corrections to the Casimir force; cf. the discussions on temperature corrections 
in [Hin].) 

The main purpose of the present note is to propose a new variant of the 
sphere-substrate configuration, namely a metal sphere suspended in a spheri- 
cally formed metallic dip. See figure 1. We will make a simple, approximate, 
calculation of the vertical Casimir force on the sphere, utilizing the known the- 
ory for the Casimir effect under conditions of perfect spherical symmetry. There 
exist several theoretical treatments of the Casimir effect under conditions of 
spherical symmetry - cf. ^1 119L I2(JI 1211 1221 for instance - but the experimen- 
tal tests of this kind of Casimir forces have so far been absent. Leaving aside 
practical difficulties such as the need of keeping the sphere in the dip in a stable 
lateral position, wc hope nevertheless that the present simple idea can be of 
interest to experimentalists. 

Temperature corrections are not included in the main formalism but are 
discussed in section 3. We use natural units, /i = c = 1, in the intermediate 
calculations, and we employ Heaviside-Lorentz units. 

2 Geometrical Set-Up 

We begin by assuming perfect spherical symmetry: let there be two concentric 
perfectly conducting singular shells situated at r = a and r — b. We shall be 
interested in the distribution of fields in the annular region a<r<b, atT~0. 
The most natural way of approach when describing this situation is to make use 
of the Green functions; for the case of spherical symmetry this kind of theory 
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Figure 1: The proposed experimental setup. A metal sphere of radius a is 
suspended in a spherically formed metallic dip of radius h. 



was worked out by Milton et al. There occur two scalar Green functions 

in the problem, Fi{r,r') and Gi{r,r'). As shown in for the double-shell 
situation, the electromagnetic boundary conditions at r = a, 6 transform into 
the following conditions for the scalar Green functions: 

FUa,r') = 0, A[^Gz(r,r')Ua = 0, (1) 

and similarly for r — b. 

The surface force density, here called on the outer surface r = 5 is calcu- 
lated from Maxwell's stress tensor. To this end we need the two-point functions, 
which are in turn given by the components of the Green functions. Some cal- 
culation yields for the radial two-point function for the electric field 



-^b^ Jo 2^ 47r s[{y) - AG{ay/b)e'i{y) 

whereas the corresponding transverse function for the magnetic field becomes 



1 , ^21 + 1 



{H^{b-)) = -— ydy 



Trb'^ .In I 471" 



si{y) - AG{ay/b)ei{y) s'i{y) - Apiay /b)e[{yy 



s'liy) ~ AG{ay/b)e[{y) si{y) - Apiay /b)ei{y) 

(3) 

Here {E^{b—)) is shorthand for {Er{r)Erir'))r'^r=b-, etc.; y means the nondi- 
mensional frequency y — ujb with uj being the Wick-rotated frequency, and 
si{x) — Iu(x), ei(x) — yj2x/i: Ki, with v = 1 + 1/2 are the Riccati- 

Bessel functions defined such that their Wronskian is W{si,ei} = — 1. Prime 
means derivative with respect to the whole argument. The coefficients Ap and 
Ag in Eq. © are 

, , , siix) , , , s',(x) 

Ap{x)^^, Ag{x)^^. 4 
ei[x) e\[x) 
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The reason why the integration over y runs to infinity in Eq. (0) is that the 

medium in the surfaces is assumed to be perfectly conducting (i.e., with per- 
mittivity e — * oo) for all frequencies. 

The other two-point functions are found to vanish, 



(i?i(M) = (ff'(M)-o, 



(5) 



so it is simple to find the surface force density at r 
stress tensor: 

f, = \U{El{h-)) + {Hl{h- 



— b via use of the Maxwell 



))] 
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1=1 



21 + 1 
47r 



s[iy)-AFiay/b)e[{y) , s['{y) ~ AGiay/b)eUy) 



si{y) ~ AF{ay/b)ei{y) s[{y) - Aciay /b)e[{y) 



(6) 

The surface force that we are interested in, is the z component of the surface 
force density integrated over the lower hemisphere, i7r<6'<7r,0<(/)< 
27r. This integration is trivial, since the magnitude of the surface force density 
contains no angular dependence. It is moreover convenient to rewrite Fz in such 
a way that the mutual contribution is separated out. Some formal manipulations 
yield 



1=1 



21 + 1 
47r 
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21 + 1 d 
4tt dy 



In 
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l-Aaix) 



4iy) 



(7) 



Here, x — ay/ b is a function of y, but when taking the partial derivative with 
respect to y in the last term, x and y are regarded as independent variables. The 
first term describes the self-force on the surface r = b, due to the fluctuating 
fields in the annular region. [If we were taking into account the contribution 
from the outer region r > b also, as would strictly speaking be necessary when 
considering a double spherical shell, then there would be an analogous term 
e'l/ei + e'l'/e'i in addition; cf. Ref. |^.] In our case, self-forces are not of 
interest. The physically important force thus becomes 



F._ 



1 

262 



y<iy 



E 

1=1 



21 + 1 d 
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In 



1-Af{x) 



eijy) 
si{y) 



1-Ag{x) 



<iy) 



(8) 

The expression is positive, corresponding to an upward directed force on the 
outer wall r = b, which in turn means a downward directed force on the sphere. 
The expression is finite as it stands; no regularization procedure is necessary. 

To make the expression more practically useful, we employ the Debye ex- 
pansion for the Riccati-Bessel functions. The calculation is parallel to that in 
Ref. |19| . and will not be repeated here. We give the result in dimensional form. 
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when reintroducing a instead of h and working to the first order in the small 
quantity d/a, 

F, = -{■Kc?)\\ + — . (9) 

Here we have separated out the standard expression 7r^?ic/240 for the Casimir 
surface force density between flat parallel plates. It is seen from Eq. that 
the curvilinear geometry leads to a slightly increased force as compared with 
the force between parallel plates having an effective area of tto^ . For instance, 
if a = 50 /im, d = \ /im, the last correction term in Eq. @ amounts to 2.6 %. 

3 Discussion 

The most characteristic property of the present proposal is that it suggests how 
the Casimir formalism worked out for spherically symmetric geometries can be 
exposed to an experimental test. As far as we know, this is the first proposal of 
such a kind. Of course, our calculation above is only approximate. Let us make 
a few final remarks: 

• The most evident simplification that we have made, is to assume that 
the field distribution is the same as in the case of perfect spherical symmetry, 
all over the dip. Of course, near = ^tt there are "stray" fields making the 
distribution different from the perfectly symmetric case. A circumstance which 
however diminishes the influence from the stray fields, is that the z component 
becomes only slightly influenced near B = ^tt. We can estimate the magnitude of 
the effect by performing the integral over 9 in the calculation from i7r + A to 
TT, instead of from ^tt to tt. The result is that the expression @ gets multiplied 
with a correction factor (1 — sin^ A). For example, taking the error introduced 
by the stray fields to be A = 5°, we get a decrease in of about 0.8 %. It 
would be quite a difficult task to make an accurate calculation of the influence 
from the stray fields; one would have to solve the complicated field distribution 
problem around 9 — ^ir. 

• It is of interest to compare our results with the new technique proposed 
by Jaffe and Scardicchio based on optical paths [231 1^ HH]' This technique 
assumes classical optics; it is most accurate at short wavelengths and where 
diffraction is not important. The technique has so far been applied to the case 
of scalar fields. The main procedure is to write the Casimir energy as a trace 
of the Green function; then the Green function is replaced by the sum over 
contributions from optical paths labelled by the number of reflections from the 
conducting surfaces. There are two central quantities in the analysis: first, there 
is the length Ir {x) of the closed geometric optics ray beginning and ending at the 
point X and refiecting r times from the surfaces; secondly, there is the so-called 
enlargement factor A^ (x) of classical optics associated with the r-reflection path 
beginning and ending at x. 

Consider first the simple sphere-plate configuration. The original wavefront 
leaving x is spherical. The first reflection from the sphere produces a new 
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wavefront, with in general two unequal radii of curvature. When next incident 
upon the sphere, the asymmetric wavefront will be transformed in a complicated 
manner, not yet worked out even for the scalar field. If we now consider our 
proposed experimental setting where there are two curved surfaces, this method 
appears to be quite complicated. We shall therefore not try to work out this, 
but it is of interest nevertheless to compare our results with those obtained for 
the sphere-plate configuration. 

Let f°P*-^'^^^{d/a) be the correction factor for a sphere and a plate, calculated 
by the optical method. This factor gives the ratio between the Casimir force 
and the force obtained for two parallel plates separated by a gap d. From 
we have, to the first order in d/a, 

joptical(^y^) = 1 + 0.05 d/a. (10) 

We can compare this with the result calculated from the proximity force ap- 
proximation (PFA) Actually there is an ambiguity in the PFA: the basic 
idea of the method is to apply the parallel-plate result to infinitesimal bits of 
the (in general) curved surfaces and integrate them up. The ambiguity lies in 
which surface is chosen for the integration. The physically best choice turns out 
to be the plate-based PFA, according to which 

/P^l°(d/a) = 1 - (l/2)d/a. (11) 

It is seen that even the signs in the correction terms in Eqs. (|10|l and Hll|l 
are different. Now, the sphere-plate situation for the scalar field has actually 
been calculated numerically, to a high accuracy [^Hl- The numerical results 
show clearly that the interaction energy increases with increasing values of d/a. 
From Fig. 4 in [21] it is seen that the optical approximation works well up to 
d/a w 0.2. 

Finally let us compare these results with our expression (9) for the Casimir 
force. We see that our expression corresponds to the correction factor 

/(d/a) = 1 + (4/3)d/a. (12) 

The experimental configuration that we have proposed in this paper is of course 
different from the sphere-plate configuration with scalar fields, but we see that 
our correction term in the factor f{d/a) is positive. There is thus a qualitative 
agreement between our result and the optical method result for scalar fields, 
Eq. (|10|l . as well as with the numerical result in |26| . 

• We next consider the correction coming from finite temperatures. This 
point is quite subtle. It is instructive to start with the case where two parallel 
plates are separated by a gap d. The general condition for applying T = 
theory with reasonable accuracy is that Td <C 1 (in natural units). Consider, 
for example, two gold plates at a gap of d = 1 /xm; from figure 5 in it follows 
that the surface pressure is about 1 mPa at T = 300 K and about 1.15 mPa 
at r = 0, thus a decrease of 15 % at room temperature as compared with zero 
temperature. As we mentioned earlier, the parallel plate experiment of Bressi 
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et al. may even have been able to measure this temperature effect for the 
case of narrow gaps, d < 0.5 /xm. [In this experiment the plates were actually 
coated with chromium rather than with gold; when d = 0.5 fiva. corresponding 
to Td = 0.065 at room temperature, the surface pressure is calculated to be 15.5 
mPa whereas the T = theory yields 20.8 mPa. The theoretically predicted 
force reduction is thus somewhat greater than above, about 25 %; cf. also the 
discussion in [TB|.] 

No consensus has so far been reached in the literature as regards the tem- 
perature correction, not even in the simple case of parallel plates. The essential 
physical point is whether the transverse electric (TE) mode contributes to the 
Casimir effect for a metal in the limit of zero frequency, corresponding to a 
Matsubara integer equal to zero. In our opinion it does not, as spelled out in 
detail in Ref. |16| . The papers of Sernelius and Bostrom are in agreement with 
this opinion |28l 1291 1801 1811 182| . It implies, as mentioned, that the Casimir 
force is weaker (by some percent when d = 1 /im) at room temperature than at 
zero temperature. By contrast, the recent paper of Chen et al. \6'6\ . which is 
based upon a reanalysis of the earlier Atomic Force Microscopy (AFM) exper- 
iment reported in ^1], claims the temperature correction to be so small as to 
be negligible. In that apparatus a gold-coated polystyrene sphere, mounted on 
a cantilever of an AFM, was brought close to a metallic surface and the deflec- 
tion of the cantilever was measured as a function of the distance. As a general 
remark on this experiment, in spite of the apparent excellent agreement of the 
experiment, one suspects that the accuracy of it has been overestimated. As 
discussed by lannuzzi et al. ) and by Milton |2] , at very short distances d (in 
the experiment the shortest distance was equal to 62 nm), the force at d = 62 
nm differs from the force at d = 62 -I- (5 by more than 3.5 pN (the experimental 
uncertainty claimed by the authors) when 5 is larger than a few angstroms. 
This implies that d should have been measured with atomic precision in order 
to correspond to the accuracy claimed. One may compare this with the real 
error in the experiment: the value of d was determined with ±1 nm accuracy. 
Moreover, the reason for the claimed smallness of the temperature correction in 
is that the analysis is based upon the plasma dispersion relation, which in 
turn implies that the zero TE mode contributes to the Casimir force. 

After these introductory remarks on the parallel-plate geometry, we now turn 
to the finite temperature version of the expression . It is convenient to work 
in terms of the free energy F, instead of the force Fz . The relationship between 
these quantities is Fz = dF/db (recall that Fz is with our sign conventions 
positive, whereas the interaction free energy F has to be negative). The finite 
temperature version of F is given by 

m=0 ;=i 

where (3 = 1/T, = / + 1/2, x = 2TTma/l3 and y = 2TTmb/(3 being the nondi- 
mensional frequencies. The prime on the m summation means that to = is to 
be taken with half weight. 



1 - ^^(a;) 



eijy) 
siiy) 



s'riy) 



(13) 
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We shall consider some limiting cases of this expression. First consider the 
static case, whereby we mean that the frequency is zero (m — 0). This means 
that we take the analytical limits of si{x) and ei{x) when a; ^ 0. The result 
becomes (cf. Ref. 22 ): 



1 °° r 
1=1 L 



(14) 



In this formula, the radii a and b are arbitrary. The case where this formula 
is most useful, is that of moderate or large gap widths. For comparison, we 
mention that at room temperature the Casimir force between a large (a = 
300 /im) gold sphere and a copper plate gets its dominant contribution from the 
zero frequency when d is greater than about 0.7 /xm (cf. Fig. 3 in '35'). The 
physical reason is that for large d it is generally the low frequencies that become 
most important. 

Consider next case of finite temperatures, assuming the gap to be narrow, 
^ = d/a <C 1. Then, from the uniform asymptotic (Debye) expansions for the 
Riccati-Bessel functions we obtain, to the first order in ^ 20 : 

oo ' oo 

/3i^ = ^ ^ J. In (l - e-2« ^^^"+^) , (15) 

m=0 1=1 

where t — 27ra//3 is the nondimensional temperature. This expression puts no 
restriction on the temperature. 

Finally let us consider the case of high temperatures. Then, only the lowest 
Matsubara frequency m — contributes, and we are led back to the expression 
H14(l . If we in addition require ^ to be small, we obtain [201 

m".^0) = ^^%d,ln(l-e-^^) = -||. (16) 

This is the same as the m = expression obtained for the free energy for parallel 
plates [33]. (Recall that when the PFA is assumed, it is this expression that is 
used for the force calculation for a sphere/plate system.) 

Numerical calculations show in general that there is for low temperatures 
a flat plateau for the free energy extending up to quite high values of t. For 
instance, when d/a = 0.05 the T = theory can be used up to logj^Q i « 1.3, or 
t = 20 (Fig. 6 in [22]). Since dimensionally t = 2TTakBT/fic, this corresponds to 
T « 140 K if a = 50 Aim. 

For low and moderate temperatures, as mentioned, the Casimir force dimin- 
ishes with increasing T. For comparison, when dealing with the force between 
a gold sphere and a copper plate, we found in [Qj the T — 300 K force to be 
weaker than the T = force by 3.6 % in the case of a width d = 0.2 //m. This is 
comparable in magnitude with the experimental results of Decca et al. [361 137) . 

To conclude: It should in principle be possible to measure the expression ^ 
for Fz, given that experimental lateral stability problems for the sphere can be 
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overcome. One should thus expect to be able to check the predicted effective 
area factor of tto^ . A practical problem is however that the last correction factor 
in Eq. ijHJl appears to be of the same order of magnitude as the temperature 
correction. 
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